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Abstract. We have carried-out a series of multi-configuration Breit-Pauli 
AUTOSTRUCTURE calculations for the dielectronic recombination of Fe^^+. We 
present a detailed comparison of the results with the high-energy resolution 
measurements reported recently from the Heidelberg Test Storage Ring by 
Schmidt et al . Many Rydberg series contribute significantly from this initial 3s^3p 
M-shell ion, resulting in a complex recombination 'spectrum'. While there is much 
close agreement between theory and experiment, differences of typically 50% in 
the summed resonance strengths over 0.1 — 10 eV result in the experimentally 
based total Maxwellian recombination rate coefficient being a factor of 1.52 — 1.38 
larger than theory over 10* - 10^ K, which is a typical temperature range of peak 
abundance for Fe^^+ in a photoionized plasma. Nevertheless, this theoretical 
recombination rate coefficient is an order of magnitude larger than that used by 
modellers to-date. This may help explain the discrepancy between the iron M- 
shell ionization balance predicted by photoionization modelling codes such as ION 
and CLOUDY and that deduced from the iron M-shell unresolved-transition-array 
absorption feature observed in the X-ray spectrum of many active galactic nuclei. 
Similar data are required for Fe*"*" through Fe^^+ to remove the question mark 
hanging over the atomic data though. 
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1. Introduction 

Dielectronic recombination (Burgess 1964) is the dominant electron-ion recombination 
process in both photoionized and coUisional plasmas. Extensive theoretical data is 
available for all K- and L-shell ions of all elements up to Zn, and selected heavy 
elements beyond, following the work of Badnell et al (2003), and is available online§. 
These data, including radiative recombination (RR), have been used to provide new 
ionization balances for both electron coUisional (Bryans et al 2006) and photoionized 
plasmas (Ferland 2006). Extensive benchmarking has taken place against experiment 
for both low-Z (C, N, O) and higher-Z elements (Fe, Ni) — see e.g. Fogle et al (2005) 
and Savin et al (2006), and references therein. Work on the M-shell is sparse (beyond 
the simple Na-like sequence, Linkemann et al 1995, Fogle et al 2003). Yet, M-shell 
Fe ions are ubiquitous in astrophysics. It has become clear recently that dielectronic 
recombination (DR) rate coefficients for Fe 3p''((; = 1 — 6) ions (Fe^+ - Fe^^"*") are 

§ Webpage |http: / /amdpp.phys.strath.ac.uk/tamoc/DATA| 
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highly questionable at temperatures where these ions form in photoionized plasmas 
(10'' — 10^ K, say, Kallman and Bautista 2001). This stems from the inability of 
photoionized plasma modelling codes such as ION (Netzer 2004) and cloudy (Kraemer 
et al 2004) to model the iron M-shell unresolved-transition-array absorption feature 
observed in the X-ray spectrum of many active galactic nuclei. The situation can be 
improved by changing the ionization balance for these Fe ions at such temperatures, 
as first suggested by Netzer et al (2003), and which is achieved by increasing the 
dielectronic recombination rate coefficients by large factors (e.g., 2 4). That this 
is a plausible approach has been verified experimentally for Fe^^+ by Schmidt et al 
(2006) who carried-out high-energy resolution DR measurements at the Heidelberg 
Test Storage Ring (TSR). They deduced a MaxwcUian recombination rate coefficient 
which is up to an order of magnitude larger than that recommended by Arnaud & 
Raymond (1992) and Mazzotta et al (1998) at photoionized plasma temperatures, 
and which is currently used to determine the ionization balance of iron in ion and 
CLOUDY, and other photoionized plasma modelling codes such as XSTAR (Kallman 
and Bautista 2001). The reason for this difference is that the existing theoretical 
dielectronic recombination contribution (Jacobs et al 1977) to the total recombination 
rate coefficient falls-off exponentially below ~ 10 eV and the total is dominated by 
direct radiative recombination. A similar problem was noted by Miiller (1999) for 
Fc"'^^"'" and it is prevclant also for L-shell ions, following the pioneering work of Savin 
et al (1997). The approach of Jacobs et al (1977) is based upon the 'no-coupling' 
scheme, allows only for dipole core-excitations in the dielectronic capture process and 
pays no detailed attention to the positioning of near-threshold resonances. It should 
be noted that the work of Jacobs et al (1977) was motivated by applications to high- 
temperature electron collision dominated plasmas and, for such, their approach is quite 
reasonable. It is clear, however, that the DR of Fe'^"^+ needs to be re-examined for 
application to photoionized plasmas. Such a re-examination, including a comparison 
with the results of the measurements by Schmidt et al (2006), will provide a benchmark 
for other Fe i'p'^iq = 1 — 6) ions, and the M-shell more generally. 

The outline of the rest of this paper is as follows: in section 2 we describe our 
theoretical approach; in section 3 we make a detailed study of the structure of Fe^^+; 
in section 4 we compare our velocity-convoluted DR cross sections with those from the 
experiment by Schmidt et al (2006); and in section 5 we compare various Maxwellian 
recombination rate coefficients. 

2. Theory 

We have used AUTOSTRUCTURE (BadncU 1987, 1997) to carry-out a series of multi- 
configuration Breit-Pauli calculations of dielectronic recombination cross sections 
and rate coefficients. The method implemented within autostructure is the 
independent processes, isolated resonances using distorted waves (IPIRDW) approach 
to DR. A detailed discussion of the validity of this approach is given by Pindzola et 
al 1992) while its advantages from a (collisional-radiative) modelling perspective is 
discussed by Badnell et al (2003). 

Let (7j^(£') denote the partial dielectronic recombination cross section, as a 
function of center-of-mass energy E, from an initial metastable state v of an ion 
through an autoionizing state j, into a resolved final state / of an ion 
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then 

a}^iE)^aj^L^{E), (1) 

where U (E) is the Lorentzian Unc shape of the resonance (energy-normahzed to 
unity). Here, a denotes the integrated (partial) dielectronic recombination cross 
section, which is given by 

where ujj is the statistical weight of the {N + l)-electron doubly-excited resonance 
state j, uJu is the statistical weight of the A''-electron target state (so, z = Z — N, 
where Z is the nuclear charge) and the autoionization (A^) and radiative {A'-') rates 
are in inverse seconds. Here, Ec is the energy of the continuum electron (with orbital 
angular momentum /), which is fixed by the position of the resonance j relative to the 
continuum v, /jj is the ionization potential energy of the hydrogen atom (both in the 
same units of energy) and (27rao)^To = 2.6741 x 10""^^ cm^s. 

A powerful aspect of the IPIRDW approach is that the use of equation 
enables an analytic integration over the resonance profiles to be carried-out. This 
is in contrast to an i?-matrix calculation which must map-out the detailed resonance 
structure numerically. This in itself is more demanding for DR than for electron- 
impact excitation since a much finer energy mesh is needed to map-out all resonances 
which contribute significantly to the cross section — see Gorczyca et al (2002) for a 
detailed study and discussion of the issue. 

So, let a^j^ denote the corresponding energy-averaged (partial) dielectronic 
recombination cross section, then 

-, nE,+AE/2 

a^iE,)^— o]^{E')dE'. (3) 

^E Je,-AE/2 ^ 

Here, /S.E denotes the bin width energy, which is chosen so as to be large compared to 
the Lorentzian width and small compared to the characteristic width of any subsequent 
convolution; otherwise, the choice of AE' is arbitrary and it is usually taken to be a 
constant. Then, 

--%m = ^^'f.m- (4) 

Thus, for a fixed j and the energy-averaged partial DR cross section takes-on a 
non-zero value at a single energy, i?c, including when summed-over final states /. 
Most applications involve a sum over resonance levels j and it is convenient to 'bin' 
the cross section via 

a,{Em)=Y.^]u{Ec) V Ec e [E„„E„+i), (5) 

where Em+i = Em + AE. The sum over / is over all final states which lie below the 
ionization limit and which may include cascade through autoionizing levels, although 
a single cascade (i.e. a two-step radiative stabilization) is usually more than sufficient. 
For total rate coefficients, applicable low-density plasmas, the sums over / and j are 
taken to convergence but for application to laboratory measurements the sum over / 
(and hence, in practice, j) is truncated. 
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2.1. Application to merged-beams measurements 



Merged-beams measurements utilizing an electron-cooler determine a rate coefBcient 
for the dielectronic recombination process. To compare with a measurement at an 
electron-ion centre-of-mass energy Eq, we determine a corresponding theoretical rate 
coefficient, a(wo), formally given by 



y.{vo) ^ <va>= / cr(w)w/(wo, v)dv , 



(6) 

where /(wo,v) is the merged-beams electron velocity distribution in the center-of- 
mass frame of the ions and vq = ^^E^jra^., since the electrons are moving non- 
relativistically with mass << mx, the mass of the ion X. 

The experimental velocity distribution, /(wo,v), is a 'flattened Maxwellian' 
(Dittner et al 1986) which is characterized by two parameters, a 'parallel' temperature 
T|| and 'perpendicular' temperature Tl, with T\\ << T±: 

(u|| - wo)^'' 



2'KkTs 



1/2 



exp 



2kTn 



27rfcTi 



■ exp 



2fcrt 



(7) 



where v\\ and denote the parallel and perpendicular components of v, respectively. 
Note, at high energies, Eq » kT±, f{vo,v) reduces to an effective Gaussian 
distribution with a full- width at half-maximum of 2(\n2 EokT^^Y^^ . 

For a bin width that is much smaller than the energy resolution of the experiment, 
and on using the distribution given by equation Q , we can write equation jnj in terms 
of the energy-averaged cross sections and bin energies, Em = 



AEa^{Em)vr, 



2kT^^ 



: exp 



So 



Er, 



X [erf {zi ~ 22) + erf {zi + Z2)] 



kTj^ — kTn 



kTj_ 



(8) 



where 



and 



Zl 



Z2 



Em{kT^~kT\\) 
kT±_kT\\ 



EpkT^ 
kT\\ (kT \ — A: Til 



1/2 



1/2 



(9) 



(10) 



Recall, a^, has units of cm^. Writing Vm/2 = y/ Em/ lay/ Ih/ (2r?T.c), we have that 
V/H/(2me) = 1.0938 x 10^ cm s ^ is the relevant remaining constant which defines 
the rate coefficient. 

In measurements carried-out at storage rings, the ions are circulated long enough, 
and the densities are low enough, for the ion population to be concentrated in the 
ground state, normally. In single-pass measurements it is necessary to calculate 
DR cross sections for metastable levels as well and then to combine them using 
experimental metastable fractions, if possible, or, typically, to use fractions for which 
the resulting cross section best matches the measured. 
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2.1.1. Survival of the species. Recombined ions with high principal quantum numbers 
are re-ionized by the strong electric field present in the charge-state analyzer which 
is used to separate the recombined ions from the original ion beam, and so are not 
counted as recombined ions. A 'hard' cut-off at ric, given by the hydrogenic expression 
(Bethe and Salpeter 1957) 

ric = (6.2 X lO^z^/i^)^^'' , (11) 

where F (V/cm) is the field strength, often suffices. Sometimes, however, recombined 
ions with n > ric have time to radiatively stabilize to n < during the time-of-flight 
(rp) from the cooler to the analyzer and so survive to be counted. This is modelled 
theoretically through the use of a 'soft' or 'delayed' cut-off (Zong et al 1998, Schippers 
et al 2001). A soft cut-off simply imposes a higher effective based-upon the lifetime 
of the Rydberg states. The delayed cut-off approach determines, for n > ric, the 
lifetime of each Rydberg state (usually hydrogenic) , Tni , and multiplies each nl partial 
DR cross section by a survival probability, given by 



„ , T„; / -Tp 

Pni = 1 cxp 



exp 1 I - exp ( 

2t„; / \ ITnl 



(12) 



where tl is the time-of-flight for the passage through the merged-beams section of the 
cooler. It is often the case that tl ^ Tni for the contributing nl and so, to a good 
approximation, 

Pra = 1 - exp : (13) 

i.e., independent of the cooler time-of-flight. The shortest lifetimes are for the lowest 
n > Uc and lowest / since the latter can radiate (n — > n') to the lowest possible n'- 
states. Thus, the flnal result is relatively insensitive to the range of n > ric considered, 
provided that there no other magnets resulting-in significant cut-offs. 

A more elaborate, experimental set-up dependent, approach is described by 
Schippers et al (2001) which takes account of the cooler geometry and the position of 
the various magnets and their fields as they impinge upon the recombined ions along 
their path to ultimate survival to be counted as recombined, or not. In addition, field 
ionization rates, due to Damburg and Kolosov (1979), are calculated explicitly. 

2.2. Application to Maxwellian plasmas 

The usual expression for the Maxwellian partial DR rate coefficient (e.g. Badnell et 
al 2003) can be obtained simply from the corresponding integrated DR cross section, 
given by equation (0): 

where (47ra^)3/2 = 6.6011 x lO'^^ cm^. Trivially, it can also be determined from the 
energy-averaged DR cross section, ct, for kT » AE, on substituting for a in equation 
p4|l from equation 

Total DR-plus-RR rate coefficients are required for plasma modelling. We 
determine the RR contribution also using AUTOSTRUCTURE, following Badnell and 
Seaton (2003) and Badnefl (2006). 
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2.2.1. Fits to totals. It is convenient often for modelling purposes to fit the total 
(Maxwellian) DR rate coefficient, a^^{T), to the following functional form: 

a?«(r) = T-3/2Vc.exp('^') , (15) 



i 

where the Ei are in the units of temperature, T, (eV or K) and the units of Ci are 
then cm'^s~^[eV or K]^/^. 

The same is also true for RR, for which we use the usual (Verner and Ferland 
1996) functional form: 

-1 



(T) = A 



(16) 



where Tq, Ti are in the units of temperature (eV or K) and the units of A are cm^s~^, 
while B is dimensionless. A more accurate representation (Gu 2003), especially for 
low-charge ions, replaces B as 

iJ^iJ + Ccxpf-^V (17) 



where, again, C is dimensionless and T2 has the units of temperature. 
3. The Fei3+ target 

The DR reactions which we take account of are defined by the iV-electron target 
configuration interaction expansion which we use. All possible {N + l)-electron 
configurations are then constructed by adding a continuum or bound orbital to them. 
All possible autoionization and (electric dipole) radiative rates are determined from 
these configurations, and are applied subsequently so as to determine partial and total 
DR cross sections, following the theory of Section |2| 

It is still the case, for Fe^'^+ at least, that it is convenient and meaningful to 
consider separately the An — and An = 1 core-excitation contributions to DR 
since, as we shall see, our highest n — i target level lies below our lowest n = 4 level. 
Indeed, 'An = 0' DR completely dominates over An = 1 at photoionized plasma 
temperatures, and can be expected to be the largest contribution too at coUisional 
plasma temperatures, at least where Fe^'^"'" is normally abundant (Mazzotta et al 
1998, Bryans et al 2006). This separation enables us to restrict the sum over Rydberg 
states to n — 1000, / = 15 and to n = 100, I — 5, for An = and An = 1, respectively. 



3.1. An = 

We consider two different target configuration interaction expansions because carrying- 
out two DR calculations enables us to contrast the level of accuracy/difference in the 
A'^-electron targets with that of differences between theory and experiment for cross 
sections and assess the accuracy of total rate coefficients, i.e. quantify the uncertainty 
in the {N + l)-electron problem. 

We define a Basis A consisting of the following configurations (assuming a closed 
Ne-like core): 

1 : Ss^Sp , 2 : SsSp^ , 3 : Ss^Sd , 4 : 3p^ 
5 : SsSpSd , 6 : Sp^Sd , 7 : SsSd^ 
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and a Basis B, which consists of Basis A plus 
8 : 3p3d2 . 

Thus, 3d'^ is the only configuration from the n = 3 complex which is omitted, by 
Basis B. Configurations 1-3, plus 5, form the minimal set which allows for all 
(An = 0) one-electron promotions during the dielectronic capture process from the 
ground configuration. Configuration 4 mixes strongly with 5, whilst 6 and 7 (which are 
strongly mixed themselves) provide the leading even parity configuration interaction. 
Configuration 8 (Basis B) provides a check on that for the odd parity. 

Basis A gives rise to 37 target terms whilst Basis B gives rise to 56 terms. In 
both cases the radial functions were determined using the Slater- Type-orbital model 
potential of Burgess et al (1989). The (3s, 3p, 3d) radial scaling parameters, A„i, were 
determined by minimizing the equally weighted sum of eigenenergies of the 18 lowest 
terms, which correspond to all of those which arise from the first 5 configurations of 
the basis expansions. For Basis A: Ass = 0.93173, Asp = 0.99255 and Asd = 0.89006. 
For Basis B: Ass = 0.94849, Asp = 1.02316 and Asd = 0.86884. All other radial scahng 
parameters were taken to be unity. 

Basis A gives rise to 84 target levels whilst Basis B gives rise to 129 levels. Our 
Breit-Pauli calculations include the one-body non-fine-structure and fine-structure 
operators, including the effective one-body Blume and Watson operator for the 
mutual-spin-orbit and spin-other-orbit interactions between valence electrons and 
the Ne-like closed shells (Badnell 1997). The effect of the two-body fine-structure 
operators representing interactions between valence electrons (including spin-spin 
now) is small — of the order 10"** Ry — and since they are time consuming to 
determine in the DR calculation we omit them, along with the two-body non-fine- 
structure operators which are of the same order effect. 

In Table we compare our lowest 40 calculated level energies obtained from 
using Bases A and B with those obtained from the NIST (2006) database and 
those calculated with superstructure by Storey et al (2000) using the Thomas- 
Fermi model potential. These levels are all of those which arise from the lowest 5 
configurations, i.e, it includes all levels which contribute to An — DR in the absence 
of configuration mixing. We note a distinct improvement in the agreement with the 
results of (Basis 2 of) Storey et al (2000), and with the observed energies, on going 
from Basis A to Basis B. Basis 2 of Storey et al (2000) included all configurations 
from the n = 3 complex as well as n = 4 configurations of the form 3s^4l and 3s3p4Z, 
for / = — 3. We note little improvement in the agreement with the observed energies 
resulting from the use by Storey et al of their larger target basis 2, compared to Basis 
B. Some high-lying levels (37-40) are now in observed order but, on the other hand, 
many of the levels of configuration 5 (3s3p3d) are distinctly higher, compared to the 
observed, than are those from Basis B. 

Nevertheless, differences of up to 0.07 Ry (mostly up to 0.03 Ry for Basis 
B) between the calculated and observed low-lying level energies means that it is 
important to use the observed target energies to position the DR resonances, so as to 
eliminate sensitivity to the exponential factor in equation H14|l at photoionized plasma 
temperatures. This is done simply by moving each (A^-t- l)-electron autoionizing level 
by the difference between the calculated and observed excitation energies between the 
initial and parent A^-electron level energies. 

Generally, A'-' << for An = DR and so a (x A'-' (see equation 12)1 . Excepting 
DR via the fine-structure core-excitation (3pi/2 ~ 3p3/2)j is dominated by the 
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Table 1. Level energies (Ry) for Fe^^+. 
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37 


5 


-2 


1 


1 


7.97268 


7.75185 


7.79103 


7.65001 


38 


5 


-2 


2 


3 


7.92298 


7.74867 


7.80163 


7.66171 


39 


5 


-2 


1 


3 


8.00254 


7.78798 


7.82900 


7.68796 


40 


5 


-2 


2 


5 


7.95856 


7.77661 


7.83024 


7.69544 



" > denotes even parity, < odd parity. 
''This work. 
Storey et al (2000). 
''NIST (2006). 



inner-electron (dipole) radiative rate. Thus, it is instructive to study radiative rates 
for Fe^^+ in some detail. In Table |2] we compare our radiative rates obtained from 
using Bases A and B with those determined by Storey et al (2000) from the 'extended' 
Basis 2 of Storey et al (1996). This 'extended' Basis 2 includes the configurations of 
Basis 2 (Storey et al 2000) but adds further configurations involving n = 4 orbitals. 
In addition, all of the n — A orbitals are now pseudo-orbitals (they were physical in 
the vanilla Basis 2) — see Storey et al (1996) for further details. As far as (total) DR 
cross sections are concerned, the distribution of radiative rates over the final states is 
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Table 2. Radiative transition 
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Basis A'' 


Basis B^ 


Basis 2*^ 


3 


1-2 


3.25(07)'= 


3.25(07) 


3.66(07) 


4 
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5 
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2.56(09) 
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3.97(10) 


13 


3-5 


4.85(08) 


1.14(09) 


1.24(09) 


13 


6-7 


3.27(09) 


3.50(09) 
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11 
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14 


10 


7.02(08) 
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^This work. 

''Extended Basis 2, Storey et al (2000). 
^=3.25(07) denotes 3.25 x lO''. 
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irrelevant, in general, so long as they are all bound. Thus, we have summed over the 
fine-structure levels of the lower term to make the comparison shown in Table |21 

Overall, we observe no drastic difference between the results obtained on using 
Basis A and B compared to those of Storey et al (2000). The results of Bases B tend 
to be closer to those of Storey et al (2000), than Basis A, for the strong radiative rates 
(~ 10^° s~^), especially from configuration 5 (3s3p3d). A general point is illustrated 
by the results for the spin-quartet level 16. This level mixes with the nearby doublet 
level 18 (both J — 3/2, odd parity). Consequently, rates from level 16 to lower- 
lying spin-doublets are very sensitive to the precise mixing. Of course, as far as DR 
is concerned, if the autoionization rates associated with these parents are such that 
A'^ » then the DR cross section is simply redistributed from one peak to another, 
the parents being less than 0.5eV apart. 

3.2. An = 1 

Art = 1 core-excitation contributions to DR come into play at high temperatures, 
i.e., in electron collision dominated plasmas. Like the case of 1 — 2 core-excitations 
in Li-, Be- and B-like ions, we expect the contribution from 'inner-shell' 2 — 3 core- 
excitations to rapidly decrease as we progress through Na-, Mg- and Al-like due to 
the increasing number of core-rearrangement autoionization channels. We consider 
both 2 — 3 and 3 — 4 An = 1 core-excitations so as to get a precise assessment of their 
relative importance and an indication of whether, or not, 2 — 3 core-excitations need to 
be considered beyond Al-like. The An = 1 contribution to a Maxwell rate coefficient 
is not sensitive to resonance positions, and so the ab initio calculated energies were 
used throughout. 

3.2.1. 3 — 4. Again, we used two different configuration basis sets for 3 — 4 core- 
excitations so as to gain insight into the uncertainty of the theoretical results. The 
first (Basis C) consisted of Basis A, plus 3s^4^ , 3s3p4^ , 3p^4^, for ^ = — 3. These 
19 configurations give rise to 250 target levels. The scaling parameters for the Slater- 
Type-Orbital model potential were determined from a subset of these configurations: 
3s^n/, for n = 3,4 and all The reason for this is to ensure that the optimization 
procedure for the n = 4 orbitals was tightly linked to the n = 4 term energies. The 
3p and 3d scaling parameters were determined first by minimizing on the lowest two 
terms. These were then fixed and the 4/ scaling parameters determined by minimizing 
the equally weighted sum of n = 4 term energies (still in the presence of the n = 3 
states). The resuh: Agp = 0.69556, Asd = 0.71607, = 0.79659, A4p = 0.7903, 
A4d = 0.7821, and A4f = 0.8563. The 3s parameter is not well determined by such a 
procedure, on the other hand, opening-up the 3s sub-shell would not be consistent 
with the procedure for the other orbitals. Hence, we simply set Ass — 0.7 for 
consistency with the other scaling parameters. Such a procedure is optimum for 
3p — > Al promotions, which can be expected to dominate the 3 — 4 core-excitations. 

The second basis that we use (Basis D) is the scattering target basis used by Storey 
et al (1996) for the electron-impact excitation of Fe^'^+. It includes all configurations 
belonging to the n = 3 complex, plus 3s^4Z , 3s3p4/ , for / = — 3. Thus, it includes the 
3p3d^ and 3d^ configurations omitted by Basis C but excludes the 3p^4Z configurations. 
The focus of their work was excitation within the n = 3 complex, including resonances 
attached to n = 4. Resonances attached to n = 4 are our primary concern here. Storey 
et al (1996) used nZ-dependent Thomas-Fermi model potentials, and the relevant 
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Table 3. Some n = 4 level energies (Ry) for Fe^''"'", relative to the ground level. 



Config. 


(25 + 1)=^ 


L 


2J 


Basis C'' 


Basis D= 


Observed'* 


4s 


2 





1 


12.9023 


13.2113 


13.0769 


4p 


-2 


1 


1 


13.9497 


14.2595 


14.2963 


4p 


-2 


1 


3 


14.0114 


14.3204 


14.3434 


4d 


2 


2 


3 


15.3886 


15.6841 


15.4549 


4d 


2 


2 


5 


15.4007 


15.6965 


15.4668 


4f 


-2 


3 


5 


16.2484 


16.5705 


16.2969 


4f 


-2 


3 


7 


16.2518 


16.5738 


16.2993 



> = denotes even parity, < odd parity. 
•^This work. 

'^Recalculated from Storey et al (1996). 
'^NIST (2006). 



Table 4. Symmetric oscillator strengths (g/) for Fe*^+ . 



Transition 


Basis C'' 


Basis D^^ 


Length 


Velocity 


Length 


Velocity 


3p - 4s 


0.2043 


0.2963 


0.3707 


0.3757 


3p - 4d 


0.9886 


1.3858 


1.5037 


1.5432 


3d - 4p 


0.1766 


0.2492 


0.2670 


0.2493 


3d - 4f 


6.0729 


7.3666 


7.4848 


7.2178 



^This work. 

'^Recalculated from Storey et al (1996). 



values of the scaling parameters are listed in their Table 1. There are 227 levels 
associated with this 17 configuration basis. 

In table|31we compare energies for the n — 4: levels which result from the dominant 
3p — > 41 promotions. We note that the 3p 4d, 4f excitation energies are better 
represented by Basis C, while Basis D is somewhat better for the lower /. In table 01 
we compare symmetric oscillator strengths for the 31 — 41' array, we use the LS-coupling 
values for simplicity. We note 20-30% differences between the results for the two bases 
— nearly a factor of 2 for 3p — 4s. There is also much closer agreement between the 
length and velocity results for Basis D compared to Basis C. Thus, it is of interest to see 
how this translates into differences in the DR cross sections. These differences affect 
not only the radiative stabilization rates but also the (dipole) dielectronic capture and 
autoionization rates, including autoionization into excited states. 

3.2.2. 2 — 3. We considered 2p — > 3^ promotions only. The contribution from 2s 31 
is < 5% of the 2p by Ne-like ions. We used a target basis which comprised the first 5 
configurations of Basis A, plus 2p^3s^3p^, 2p^3s^3p3d, and 2p^3s3p'^; the latter so as 
to allow for the strong configuration mixing with the prior. These configurations give 
rise to 178 target levels. In addition, we now need to allow for core-rearrangement 
autoionization transitions of the form: 2p^3l^nl' — > 2p^3Pnl' + e~, i.e., where the 
Rydberg electron is a spectator. These transitions strongly suppress DR since they are 
'additional' autoionization pathways (ones which have no reverse dielectronic capture 
process to balance them) and are independent of n, while the populating dielectronic 
capture rate scales as n~^. 
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4. Velocity-convoluted DR cross section results 

In order to make a comparison between our theoretical DR cross sections and those 
measured by Schmidt et al (2006) at TSR, we convolute them according to equation 
© using the experimentally determined kT±^ = 1.2 x 10^^ eV and fcT|| = 9.0 x 10^^ eV, 
and apply initially a hard cut-off (see eauation lll|l at nc = 45. 1| It is assumed that there 
is no significant metastable fraction in the ion beam when it comes to the comparison 
with experiment, and so we consider only DR from the ground level of Fe^'^"''. 

Parental contributions 

Since many Rydberg series contribute to the final observed 'spectrum', we first show 
results for various core-excitations, labelled according to parent level (see table or 
configuration. We also compare results obtained using target Bases A and B. Note 
the caveat, parentage is not a good quantum number. 

In Figure ^ we show the convoluted theoretical DR results, which have the 
dimension of a rate coefficient, for the fine-structure core-excitation, i.e., parent level 

I to parent level 2, as listed in Tabled This opens- up at n = 32 and falls-off rapidly 
in n, as radiative stabilization takes place via an outer-electron transition n n' , for 
n' < 32, and which is dominated by n' < 10. The results from Basis A and Basis B 
are indistinguishable in the figure. 

In Figure |21 we illustrate the complexity which arises from the 3s — > 3p core- 
excitation to parent configuration 2. The lowest spectrum is the sum of contributions 

II The value of nc = 55 given by Schmidt et al (2006) is incorrect (Schmidt, private communication). 
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Figure 2. Velocity-convoluted DR cross sections for the 3s — > 3p core-excitation 
to parent configuration 2. Bottom curve, summed-over all parent levels, i = 3 — 10. 
Offset by 1 X 10~^ -I- (i — 3)5 X 10"^", the contributions from the individual parent 
levels, i. 



from the ^Pj, ^Dj, ^S,j and ^Pj levels. Offset above are the individual parent level, 
i, contributions for i = 3 — 10. The results from Basis A and Basis B are barely 
distinguishable in the total spectrum. 

In Figure |31 we present results for the 3p —^ 3d core-excitation. There are 
only two parent levels and the ^05/2 parent level 12 contributes only weakly as the 
J = 5/2 1/2 core radiative transition is electric dipole forbidden. Again, the results 
from Basis A and Basis B are indistinguishable in the figure. 

In figure^ we compare the results from Basis A and Basis B for excitation of the 
parent 3p'^ configuration 4, summed-over all 5 parent levels. This excitation only takes 
place through configuration mixing and we note (i) that it is weaker, especially above 
10 eV, and (ii) that we now see some small differences between the results of Basis A 
and Basis B. The lowest autoionizing states have n = 6. Core radiative stabilization 
is allowed to parent configuration 2 via 3p 3s. These 'recombined' levels then first 
autoionize at n = 7 or 8, for the spin-doublet parents. 

In figure |5l we compare the results from Basis A and Basis B for the 3s 3d 
core-excitation of parent configuration 5, summed-over all 23 parent levels. Again, 
there are small differences, below 15 eV. The lowest autoionizing states have n = 5 or 
6, depending on the parent. There are two main core radiative stabilization pathways: 
3d — > 3p, to parent configuration 2, and 3p — » 3s, parent configuration 3. These 
then first autoionize at n = 7 for the latter, and, again, between n — 7 and 8 for 
spin-doublet parents of the former. So, just a few n-values contribute strongly, but 
they are spread out in energy because of the 25 eV spread of levels of configuration 5. 

In figures El to |H1 we present results for core-excitations to parent configurations 
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Figure 3. Velocity-convoluted DR cross sections for the 3p — > 3d core- 
excitation to parent configuration 3. Bottom curve, summed-over all parent levels, 
i = 11 — 12. Offset by (i — 10)2 X 10~^, the contributions from the individual 
parent levels, i. 
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Figure 4. Velocity-convoluted DR cross sections for excitation of the parent 3p'^ 
configuration. Solid (red) curve, Basis A; dashed (green) curve. Basis B. 
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Figure 5. Velocity-convoluted DR cross sections for the 3s — > 3d core-excitation 
of parent configuration 5. Solid (red) curve, Basis A; dashed (green) curve, Basis 
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Figure 6. Velocity-convoluted DR cross sections for excitation of the parent 
configuration 6. Solid (red) curve. Basis A; dashed (green) curve. Basis B. 
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Figure 7. Velocity-convoluted DR cross sections for excitation of the parent 
configuration 7. Solid (red) curve, Basis A; dashed (green) curve, Basis B. 
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Figure 8. Velocity-convoluted DR cross sections for excitation of the parent 
configuration 8. Solid curve, Basis B. 
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Figure 9. Velocity-convoluted recombination cross sections for Fe^^+ near 
threshold. Solid (red) curve, experimental results of Schmidt et al (2006); long- 
dashed (green) curve, theoretical DR results obtained on using Basis A; short- 
dashed (blue) curve, theoretical results obtained on using Basis B; dotted (purple) 
curve, theoretical RR results (Basis A). Theoretical results are all this work. 

6-8, which are only accessible via configuration mixing. Again, there is an allowed 
core-radiative stabilization pathway for limited n-values. We see that the DR cross 
sections are progressively weaker as we go to higher-energy parent configurations. The 
main difference now between the results of Basis A and Basis B, see figures El and [7| is 
an energy shift, due to the fact that we only adjusted to the lowest 40 observed level 
energies. 

We conclude that the use of observed energies mitigates against differences in 
the DR cross section due to the different level energies of Basis A and Basis B while 
differences in the radiative rates either occur for transitions which do not contribute 
strongly to the DR, or the rates themselves are simple redistributed amongst near-by 
levels. 

4-2. Comparison with experiment 

We consider the An — and An = 1 core-excitations separately. (Low-lying 
resonances which arise from n = 3 — 4 capture to n = 4 do overlap the An — 
energy range but the peaks are so small as to be 'lost in the noise' when comparing 
with experiment.) 

4-2.1. An = 0. In figure we compare our theoretical DR results, obtained using 
Basis A and Basis B, with the experimental measurements of Schmidt et al (2006). 
The comparison is made very close to threshold (10^"* — 1 eV), utilizing a log-log scale. 
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Figure 10. Velocity-convoluted An = recombination cross sections for 
Fe^^+ over — 2.5 eV. Solid (red) curve, experimental results of Schmidt et 
al (2006); long-dashed (green) curve, theoretical DR results obtained on using 
Basis A; short-dashed (blue) curve, theoretical results obtained on using Basis B. 
Theoretical results are all this work. 

Schmidt et al (2006) note a relatively small anomalous enhancement below 10"'^ eV 
and so estimate the 'true' DR-plus-RR contribution to tend towards 1 x 10~^ cm'^ s~^ 
at 10~^ eV. The two sets of theoretical results are in accordance above ^ 0.01 eV 
while at lower energies the results of Basis B are in somewhat better agreement with 
experiment, down to ~ 0.001 eV. Below the dip at about 0.04 eV, Schmidt et al show 
fits to 5 DR resonances. We find 14 resonances, of the form 3s^3p3/232Z(^ — 0, 1), 
3s3p^*P3/29f, and 3s3p3d5g. Of course, our convoluted cross sections only exhibit 3 
(Basis A) or 4 (Basis B) obvious peaks in this energy region. We have looked at the 
bound states just below threshold but they are well below the apparent uncertainties 
in resonance energy positioning seen in figure El There is substantial disagreement 
between theory and experiment between about 0.03 and 0.35 eV. Resonances in this 
region are higher members of those just series mentioned, as well as those attached to 
parent configurations 6 and 7. However, the differences below 0.2 eV have negligible 
effect on any differences in the Maxwellian rate coefficient at 2 eV. Only the differences 
up to 0.35 eV start to impact upon the Maxwellian rate coefficient at 2 eV. We also 
illustrate the RR contribution in this energy region, having applied a hard cut-off at 
ric = 45 again. 

In figure 1101 we compare the results of theory and experiment near threshold 
again, this time using a linear plot. We see more clearly now that the theoretical 
results from Basis A and Basis B are in close agreement over « 0.05 — 0.4 eV, but differ 
substantially from the measured. Overall, in this energy range, the differences between 
the two sets of theoretical results are not large enough to suggest an uncertainty 
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Figure 11. As figure 10, but over 2 — 10.5 eV. 

which could account for the difference with the measured. Although, where there 
are more noticeable differences, the results from Basis B are perhaps somewhat of an 
improvement over Basis A. The DR cross section drops substantially (by a factor of 
10, or so) above 2.34 eV because autoionization into the excited fine-structure level 2 
becomes energetically allowed. This means that resonances below 2.34 eV contribute 
'disproportionately' at higher temperatures — see, e.g., figure 4 of Schmidt et al 
(2006). If we sum-up the resonance strengths over 0.35 to 2.5 eV we find that the 
result for Basis A is 3% smaller than for Basis B but the measured is 33% larger. 

In figure 1111 we compare the results of theory and experiment over 2 — 10 eV. 
Although this is the typical temperature range of peak abundance for Fe^^+ in a 
photoionized plasma, only above a temperature corresponding to about 5 eV is the 
Maxwellian rate coefficient dominated by the resonances above 2 eV. Again, the 
differences between the results of Basis A and Basis B are not too significant, but the 
agreement with the measured is rather poor. If we sum-up the resonance strengths 
over 2.5 to 10.5 eV we find that the result for Basis A is 8% smaller than for Basis B 
but the measured is 40% larger. Thus, over 2 — 10 eV we expect the experimentally 
deduced Maxwellian rate coefficient to be roughly 40% larger than the theoretical one. 

In figure 1121 we compare results over 10 — 20 eV. As we move up in energy we 
start to see more of a convergence between the theoretical results and the measured; 
the summed resonance strengths from Basis A and Basis B differ by less than 1% 
whilst the measured is 21% larger. However, if we now look at the energy range 
20 — 30 eV (figure [T^ we see that the measured rate coefficient appears to be sitting 
on a much larger 'background' over 23 — 29 eV — the theoretical results drop much 
closer to zero between the main resonance peaks, i.e., there is little 'fill-in' due to 
other small resonance contributions. Over 30 — 40 eV (figure I14|) the much better 
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Figure 12. As figure [1771 but over 10 - 20 eV. 
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Figure 13. As figure IT?1 but over 20 - 30 eV. 
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agreement between theory and experiment is resumed, except over 36 — 37 eV. Here, 
there is a noticeable contribution from capture to high-n states and the comparison 
with experiment is dependent on (the modcUing of) their survival to be detected. We 
illustrate with the results fr'om two models (both for Basis A) and discuss them in 
detail next. 

In figureslTslandllGlwe compare results over 40—50 eV and 50—60 eV, respectively. 
The resonances arise from 3s —>■ 3p and 3p — > 3d core-excitations for the former, but 
the second only for the latter — see figures [21 and |31 The results of Basis A and 
Basis B are indistinguishable on this scale, and so only Basis A results are shown 
and considered further. In these two energy regions we address the issue of the 
survival of the recombined states as the ions travel from the cooler to the charge- 
state-analyzer. Recall figure 121 there are high-n DR contributions from parent levels 
8 — 10 of configuration 2 which span 44 — 49 eV. (Those attached to parent level 6 
affect the 36 — 37 eV range.) The situation, as illustrated by figure O is simpler for 
the peak at 58 eV. It is in these two energy ranges for which a hard cut-off at ric — 45 
results- in large discrepancies between theory and experiment. The discrepancy is 
reduced on implementing a delayed cut-off (equation 1131) utilizing the appropriate 
time-of-flight for this experiment of 166.5 ns and imposing a hard cut-off at Uc = 95 
due to the correction magnet close to the cooler (Schippers, private communication). 
Further improvement in agreement is obtained on utilizing the survival probabilities 
of Schippers et al (2001), but for the DR of Fc"'^'^"'' . These latter two sets of results are 
only shown at energies where they differ from the 'hard cut-off' results. Clearly, the 
final agreement between theory and experiment is sensitive to the precise contribution. 
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Figure 17. Velocity-convoluted DR cross sections for 3 — 4 core-excitations in 
Fe^^+ . Bottom curve, summed-over all parent configurations. Offset (increasing) 
are the contributions from 3p — > 4s, 4d; 3s — > 4p, 4f; 3p — > 4p, 4f ; and 3s — > 4s, 4d 
core-excitations . 



i.e. survival, of recombined states with n > 45. While this largely accounts for the 
discrepancies over 44 — 49 eV between experiment and theory utilizing only a hard 
cut-off, as well as the mis-match in the final position of the Rydberg peak at 58.5 eV, 
it has no effect on the puzzling discrepancy between about 56 and 58eV. Here, the 
experimental result actually lies at a fairly uniform 4 x 10^^" cm'^s"^ below (all of) 
the theoretical one(s). 

We close the discussion of the An — results with a small observation: the 
noticeable drop in the DR cross section just below 56.5 eV (in figure IT^ is due to the 
final-state of the 3d3/2 — > 3p3/2 core-radiative stabilization pathway infact opening-up 
at this point (n = 32) to autoionization to the 3pi/2 continuum. (We obtain 0.220 
for the ratio of the 3d3/2 — >■ 3p3/2 to 3d3/2 — > 3pi/2 radiative rates compared to 0.225 
obtained by Storey et al (2000) — see also table |21) Thus, there appears to be close 
qualitative agreement between theory and experiment for this effect. This is in stark 
contrast to the poor quantitative agreement for the absolute cross section. 



4-2.2. An — 1. In figure [T7I we present an overview of the different contributions to 
the complete 3 — 4 'spectrum'. The first point to note is that the resonance strengths 
are now a factor 10, or more, smaller than those we have seen associated with the 
An = core-excitations. Only the 3p 4s, 4d core-excitations exhibit the classic DR 
spectrum. The dipole core-excitations dominate, along with 3s 4d. We see that the 
final total spectrum is complex, which makes it difficult to identify individual peaks 
in the measured spectrum. 

In figures [THll21l we make comparisons between the results of our calculations with 
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Figure 18. Velocity-convoluted An = 1 (3 — 4) DR cross sections for Fe^^+ 
over 60 — 100 eV. Solid (red) curve, experimental results of Schmidt et al (2006); 
long-dashed (green) curve, theoretical DR results obtained on using from Basis 
C; short-dashed (blue) curve; theoretical results from Basis D. Theoretical results 
are all this work. 



those of the experimental measurements by Schmidt et al (2006) for DR in the energy 
region in which the 3 — 4 core-excitations contribute. 

In figure 1181 it is capture to n = 4 which dominates, and we did not make a 
separate calculation using Basis D since the orbitals were based on the Thomas-Fermi 
potential optimized for the Al-like core. In the case of Basis C, although we used the 
same scaling parameters for capture to n = 4 as for n > 4, the Slater- Type-Orbital 
model potential depends on configuration (actually, the complex) and so is inherently 
'optimized' differently for capture to n = 4 and n > A (see Burgess et al 1989 for 
specific details). Hence, we added the Basis C n — A results to those of Basis D. In 
the 60 — 90 eV range the theoretical cross sections are much weaker than the measured. 

In contrast, in figures [T^ and 1^ we see that the theoretical DR cross sections are 
more strongly peaked than the measured, especially so for Basis D which has the larger 
4 — > 3 radiative rates. Despite Basis C apparently giving rise to a worse structure 
for Fe^'^"'" than Basis D, based upon the agreement between the length and velocity 
forms for the relevant oscillator strengths (see table EJ, it does appear that Basis C 
gives rise to distinctly better agreement with experiment for the DR cross sections 
than Basis D. Although the agreement is worse than for An = core-excitations, the 
sensitivity of the atomic structure to the overlap of the n — i and n — A orbitals 
means that the differences are 'less serious', i.e., there is still enough uncertainty in 
the atomic structure so that the agreement might be improved upon. We note that 
there is no simple identification of the various peaks in these figures as many different 
core-excitations contribute in this energy region, but they are dominated by a few low- 
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Figure 19. As figure [TSl but over 100 - 140 eV. 
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Figure 20. As figure ITsl but over 140 - 180 eV. 
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Figure 21. As figure [TSl but over 180 - 230 eV. 



n resonances because An — autoionization of the final states opens-up at n < 10 
for all parents, except 3s^3p. In addition, the 3p 4s series Rydberg accumulation 
contributes to the peak at 175 eV. 

In figure 1^ we can unambiguously identify the peaks around 210 and 220 eV as 
being associated with the limit of the Bp — > 4d and 3s 4p core-excitations, but the 
theoretical cross sections are much smaller than measured for these peaks, especially 
so for the latter. The sum over n is fairly well converged by n = 45 and even applying 
no cut-off (or infinite time-of-flight) does not increase the size of the theoretical peaks 
by much. We do see the 'overshoot' of the Basis D results here — as expected from 
the use of the unadjusted energies of Table |31 

Summing over all DR resonance strengths for the 3 — 4 core-excitations we find 
the Basis D result to be 50% larger than that for Basis C, which is comparable with 
the excess of the 3p — (4s + 4d) oscillator strength seen in Table ^ The sum of the 
measured DR resonance strengths in the 60 — 240 eV range is 21% larger than that from 
Basis C. Of course, we have noted significant disagreements between the calculated 
and measured DR resonances, both over- and under-estimates. 

Finally, in figure we present our results for the 2 — 3 core-excitations. These 
resonances lie above the highest energy considered by Schmidt et al (2006). It should 
be noted, however, that the sum of DR resonance strengths associated with this 2 — 3 
core-excitation is a factor of 2.5 larger than that associated with the 3 — 4. Thus, 
apart from contributing at a higher temperature, this core-excitation is likely to be 
more important than the 3 — 4 for application to coUisional plasmas. 
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Figure 22. Velocity-convoluted An = 1 (2 — 3) recombination cross sections for 
Fe^^+. Solid curve, present theoretical results. 



5. Maxwellian rate coefficients 

In figure we present our theoretical results for Maxwellian rate coefficients: RR- 
plus-DR from 3 — 3, 3 — 4 and 2 — 3 core-excitations, and compare the sum total of 
these with the one determined by Schmidt et al (2006), based primarily upon their 
measured DR (cooler) rate coefhcients. Over lO'* — 10^ K, a typical temperature range 
for photoionized plasmas where Fe^'^"'" is abundant, the experimentally based total is 
between a factor 1.52 — 1.38 larger than our theoretical total — this is inline with what 
we expect following our earlier detailed comparison of the DR resonances contributing 
at these temperatures. Nevertheless, it is clear that the total recombination rate 
coefficient of Fe^^+ in photoionized plasmas is an order of magnitude larger than has 
been used to-date, as first pointed-out by Schmidt et al (2006) on the basis of their 
measurements for this ion. We show also only the low temperature fall-off of the 
recommended DR rate coefficient of Arnaud and Raymond (1992) and not any of the 
ad hoc changes proposed by Netzer (2004) and by Kraemer et al (2004). 

At temperatures of a few times 10^ K, typical of electron collision dominated 
plasmas where Fe^^+ is abundant, the experimentally based total is only 5% smaller 
than our calculated one, while the recommended data of Arnaud and Raymond (1992) 
lies only about 10% higher. (The data of Arnaud and Raymond (1992) is based 
principally upon the results of Jacobs et al (1977), but includes an estimate of the 
contribution from 2p— 3d inner-shell transitions as well and which were not included by 
Jacobs et al .) We see also that both An = 1 contributions contribute only modestly 
to the total, both equally about 6% at 3 x 10^ K. Far off equilibrium, this rises to 
about 20% at 10^ K with three-quarters coming from the 2 — 3 core-excitation. Given 
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Figure 23. MaxwcUian rate coefficients for Fe^''+. Solid (red) curve, total DR- 
plus-RR; short-dashed (blue) curve, RR; long-dashed (green) curves, DR for 3 — 3, 
3 — 4 and 2 — 3 core-excitations. All this work. Dot-dashed (light blue) curve, DR 
of Arnaud and Raymond (1992). Dotted (purple) curve, experimentally based 
total of Schmidt et al (2006). PP and CP denote typical photoionized and 
electron coUisional plasma temperature ranges, respectively, for Fe^^+ (Kallman 
and Bautista 2001 and Mazzotta et al 1998). 



the modest contribution from An = 1 core-excitations, then, since such calculations 
are more demanding than for An = 0, the use of LS'-coupling may suffice, even with 
the 30% difference from intermediate coupling which we note (not shown). 

These rate coefficients are for the 3s^3p^Pj=i/2 ground level of Fe^'^+. In 
photoionized plasmas, the ion population may not be concentrated in the ground 
level, while in electron coUisional plasmas there may be signification population of 
levels of the 3s3p^ "^P term. The rate coefficients for the various metastable levels can 
be quite different. In figure !^ we compare and contrast total DR (i.e., summed-over 
all core-excitations) and RR rate coefficients from the ground and metastable levels. 
We index the target metastable levels by m, where rn = 1 for the (J = 1/2) ground 
level, rn = 2 for the other (J = 3/2) fine-structure level and rn = 3,4,5 for the 
(J — 1/2, 3/2, 5/2) levels of the metastable term, as per tabled 

At photoionized plasma temperatures we see that the m = 2 metastable DR rate 
coefficient is an order of magnitude smaller than for the (m — 1) ground level — this is 
due primarily, of course, to the absence of the fine structure DR pathway. The DR rate 
coefficients for higher metastables exhibit irregular behaviour (at low temperatures) 
due to the positioning close to threshold of the lowest autoionizing states, relative to 
these excited Fe^'^+ thresholds, although the rn = 3 and 4 metastable levels again have 
fine-structure dielectronic capture pathways and the enhancement for to = 3 appears 
to be quite noticeable. At coUisional plasma temperatures the DR rate coefficients 



Dielectronic recombination of Fe^^^ 



29 




Table 5. DR fitting coefficients (cm'^s ^K'^/-^) and Ei{K) for tfie ground and 
mctastable levels (m = 1 — 5) of Fe^^+ . 



m 


Cl 


C2 


C3 


C4 


C5 


C6 


C7 


C8 


1 


1.090(-3) 


7.801(-3) 


1.132(-2) 


4.740(-2) 


1.990(-1) 


3.379(-2) 


1.140(-1) 


1.250(-1) 


2 


3.176(-4) 


1.097(-3) 


1.451(-2) 


4.623(-2) 


1.424(-1) 


3.105(-2) 


1.173(-1) 


1.579(-1) 


3 


9.230(-4) 


4.787(-3) 


7.598(-3) 


1.538(-2) 


1.512(-2) 


1.711(-2) 


9.083(-3) 


4.875(-l) 


4 


6.837(-4) 


3.386(-3) 


8.737(-3) 


2.334(-2) 


2.819(-2) 


1.282(-2) 


9.735(-3) 


2.670(-l) 


5 


5.606(-4) 


3.306(-3) 


9.372(-3) 


1.635(-2) 


1.674(-2) 


1.783(-2) 


9.195(-3) 


3.689(-l) 


m 


El 


E2 


E3 


E4, 


E5 


Ee 


E7 


Es 


1 


1.246(3) 


1.063(4) 


4.719(4) 


1.952(5) 


5.637(5) 


2.248(6) 


7.202(6) 


3.999(9) 


2 


1.204(3) 


1.214(4) 


5.689(4) 


1.983(5) 


5.340(5) 


2.414(6) 


7.302(6) 


4.245(9) 


3 


9.887(2) 


1.075(4) 


4.542(4) 


1.832(5) 


5.506(5) 


1.692(6) 


6.994(6) 


2.789(9) 


4 


1.824(3) 


1.101(4) 


4.775(4) 


2.131(5) 


1.109(6) 


5.479(6) 


1.365(9) 


2.517(9) 


5 


1.288(3) 


1.171(4) 


4.780(4) 


1.828(5) 


5.530(5) 


1.697(6) 


6.960(6) 


2.666(9) 



split primarily into two groups which are based upon term, not level, as the influence 
of fine-structure DR and threshold effects is diminished. 

The difference in RR rate coefficients is much less pronounced, at all temperatures. 
At low temperatures the stability against autoionization for all n distinguishes RR of 
the ground level. There is little difference for metastables levels of the excited term 
because the highest stable recombined n is (almost) independent of the fine-structure 
parent. 

In tables |S1 and we present separately the fitting coefficients for our total DR 
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Table 6. RR fitting coeflicients for the ground and metastable levels (m 
of Fei3+. 



m 


A (cm^E 




B 


To(K) 


Ti(K) 


C 


T2(K) 


1 


4.321{- 


-10) 


0.6091 


2.255(03) 


4.962(07) 


0.0356 


1.006(05) 


2 


2.031{- 


-11) 


0.5464 


2.669(05) 


5.310(07) 


0.0277 


9.907(08) 


3 


1.591(- 


-09) 


1.0274 


1.196(01) 


2.038(07) 


0.0449 


1.764(08) 


4 


1.591{- 


-09) 


1.0274 


1.196(01) 


2.038(07) 


0.0449 


1.764(08) 


5 


4.803{- 


-11) 


0.3781 


1.920(04) 


4.025(07) 


0.5220 


5.410(03) 



and RR rate cocfBcients for the ground and metastable levels (indexed by m), which 
are based upon the functional forms given by equations (|15|1 and 116|l . respectively. 
The fits are accurate to better than 1% over z^(10^ — fO^) K, where z — 12> here. 



6. Summary 

We have carried-out a series of multi-configuration Breit-Pauli calculations for the 
dielectronic recombination of Fc"'^'^"''. Whilst there is much agreement between the 
theoretical velocity-convoluted cross sections and those determined experimentally by 
Schmidt et al (2006), differences over 0.1 — 10 eV lead to the experimentally based 
total Maxwellian recombination rate coefficient being upwards of 50% larger than the 
theoretical one over the temperature range 10** — 10^ K, which is typical of photoionized 
plasmas where Fe^^+ is abundant. Such a difference lies well outside of the theoretical 
uncertainty, based-upon the accuracy of the radiative rates and sensitivity to resonance 
positions. It is also well outside of the experimental uncertainty of ±18% (Schmidt et 
al 2006). 

It is difficult to see how to resolve this difference. Simply carrying-out a larger 
(configuration interaction) calculation would not be expected to result in a change 
much beyond the difference, already noted, between the Basis A and Basis B results, 
especially given the level of agreement for radiative rates which we have observed 
between Bases A and B and the extended Basis 2 of Storey et al (1996). Furthermore, 
the good agreement between theory and experiment at higher An = energies is 
interrupted twice by disconcerting differences spanning several electron volts. Such 
differences are also outside of the range of higher-order effects such as interacting 
resonances and the interference between DR and RR. Perhaps the results of a separate, 
independent, calculation will shed some light on the matter. 

Nevertheless, the theoretical recombination rate coefficient determined here for 
Fe^'^+ is an order of magnitude larger than has been used by modellers in the past. 
This may help explain the discrepancy between the iron M-shell ionization balance 
predicted by photoionization modelling codes and that deduced from the iron M- 
shell unresolved-transition-array absorption feature observed in the X-ray spectra of 
many active galactic nuclei. New data are clearly required for the other Fe Sp"? ions, 
especially g = 2 — 5, in order to eliminate the uncertainty in the DR atomic data used 
by CLOUDY, ION and xstar and to enable them to focus on the 'bigger picture'. 
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